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ABSTRACT 

A  new  ocean  mixed  layer  model  (OMLM)  was  embedded  into  an  ocean  general  circulation  model  (OGCM) 
with  the  aim  of  providing  an  OGCM  that  is  ideal  for  application  to  a  climate  model  by  predicting  the  sea  surface 
temperature  (SST)  more  accurately.  The  results  from  the  new  OMLM  showed  a  significant  improvement  in  the 
prediction  of  SST  compared  to  the  cases  of  constant  vertical  mixing  and  the  vertical  mixing  scheme  by  Pacanowski 
and  Philander.  More  accurate  prediction  of  the  SST  from  the  new  OMLM  reduces  the  magnitude  of  the  restoring 
term  in  the  surface  heat  flux  and  thus  provides  a  simulated  ocean  that  can  be  coupled  to  the  atmospheric  general 
circulation  model  more  naturally.  The  new  OMLM  was  also  shown  to  improve  various  other  features  of  the 
OGCM  such  as  the  mixed  layer  depth  and  the  equatorial  circulation. 


1.  Introduction 

Extensive  efforts  have  been  devoted  to  the  advance 
of  the  ocean  general  circulation  model  (OGCM)  during 
last  few  decades,  with  the  intention  of  simulating  the 
present  climate  and  predicting  the  future  climate.  For 
the  purpose  of  climate  prediction,  the  OGCM  must  be 
coupled  to  an  atmospheric  general  circulation  model 
(AGCM),  along  with  the  properly  represented  air-sea 
interaction.  This  particularly  requires  accurate  predic¬ 
tion  of  the  sea  surface  temperature  (SST)  from  the 
OGCM.  Hence,  we  should  regard  the  accurate  predic¬ 
tion  of  SST  as  the  most  important  prerequisite  of  the 
OGCM  for  its  application  to  climate  prediction. 

From  the  early  stage  of  OGCMs  the  surface  boundary 
condition  of  temperature  in  OGCMs  has  been  given 
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predominantly  by  the  restoring  boundary  condition  in 
which  the  simulated  SST  is  restored  toward  the  cli¬ 
matological  SST.  That  is,  the  surface  heat  flux  Q  is 
given  by 

Q  =  pC,Azm  -  Ts)!t,  (1) 

or  it  is  incorporated  into  the  equation  of  temperature  in 
the  first  layer  as 


Here  t  is  the  restoring  timescale,  Ts  is  the  temperature 
of  the  first  model  layer,  Tf  is  the  observed  SST,  p  is 
the  density,  Cp  is  the  specific  heat  of  seawater,  and  Azj 
is  the  thickness  of  the  first  layer  of  the  model. 

The  restoring  boundary  condition  has  been  favored 
in  particular  since  it  insures  that  the  simulated  SST  will 
not  drift  far  from  the  climatological  SST  and,  thus,  helps 
to  reproduce  an  ocean  that  is  close  to  the  climatological 
ocean,  even  in  the  presence  of  considerable  defective¬ 
ness  of  the  model.  Its  widespread  use  is  also  attributed 
to  the  fact  that  there  were  no  available  reliable  clima- 
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tological  heat  flux  data  and  that  the  simulation  solely 
based  on  the  prescribed  heat  flux  suffers  from  drift  from 
the  climatological  values  with  time  because  of  the  ab¬ 
sence  of  feedback  (Cane  1994;  Seager  et  al.  1995; 
McWilliams  1996). 

Nonetheless,  the  restoring  boundary  condition  im¬ 
poses  various  serious  problems  for  realistic  simulation 
of  the  ocean  using  an  OGCM.  The  most  conspicuous 
defect  is  that  the  accurate  prediction  of  SST  leads  to 
the  unrealistic  situation  of  no  heat  flux  at  the  sea  surface. 
It  also  diminishes  the  amplitude  of  the  annual  cycle  and 
causes  a  phase  lag  compared  to  the  climatological  SST 
(Haidvogel  and  Bryan  1992).  Meanwhile,  it  leads  to 
weakened  convective  overturning  and  thus  reduced  ther¬ 
mohaline  circulation  since  it  tries  to  maintain  the  cli¬ 
matological  SST  (Mikolajewicz  and  Maier-Reimer 
1994;  Cai  and  Chu  1996;  Toggweiler  et  al.  1989).  Ac¬ 
cordingly,  it  was  pointed  out  that  the  so-called  ther¬ 
mohaline  catastrophe  may  be  related  to  inadequacy  of 
the  restoring  boundary  condition  (Cai  and  Godfrey 
1995;  Rahmstorf  and  Willebrand  1995).  Moreover, 
strong  artificial  damping  to  the  climatological  SST  sup¬ 
presses  the  variability  at  the  sea  surface,  which  must 
originate  naturally  from  the  internal  ocean  dynamics, 
and  thus  leads  to  the  serious  underestimation  of  sea 
surface  height  variability  and  eddy  kinetic  energy  in  the 
high-resolution  OGCM  compared  to  the  satellite  data 
(Stammer  et  al.  1996;  Ishikawa  et  al.  1997). 

The  restoring  boundary  condition  of  the  OGCM  is 
particularly  inappropriate  in  a  climate  model  in  which 
heat  flux  is  calculated  directly  from  the  AGCM  instead 
of  by  the  restoring  term.  Relaxation  of  the  restoring  term 
for  the  coupling  to  the  AGCM  leads  to  the  generation 
of  spurious  oscillatory  behavior  in  the  OGCM  (Cai  and 
Chu  1996).  It  may  be  caused  by  the  fact  that  the  re¬ 
storing  boundary  condition  forces  the  simulated  result 
toward  an  unnatural  condition  dictated  by  the  clima¬ 
tological  SST,  thus  making  it  inconsistent  with  the  in¬ 
ternal  ocean  dynamics  controlled  by  the  model. 

Another  problem  regarding  the  restoring  boundary 
condition  employed  in  OGCMs  involves  the  magnitude 
of  the  restoring  timescale  r.  The  restoring  timescale  used 
in  present  OGCMs  is  typically  about  30  days  (Semtner 
and  Chervin  1992;  Stammer  et  al.  1996;  Marotzke  and 
Willebrand  1991).  However,  analyses  of  the  observed 
data  revealed  that  sensitivity  of  the  surface  heat  flux  to 
variations  of  SST  must  be  much  weaker  than  presumed 
in  the  above  values,  which  implies  that  t  must  be  much 
longer  (Seager  et  al.  1995;  Barnier  1998;  Chu  et  al. 
1998).  For  example,  Barnier  (1998)  and  Frankigonoul 
et  al.  (1998)  suggested  that  the  restoring  timescale  must 
be  at  least  twice  as  long. 

Recently,  improvement  of  the  restoring  boundary 
condition  has  been  attempted  by  combining  the  restoring 
term  with  the  climatological  heat  flux  Q*.  That  is, 

Q  =  Q*  +  pCrAz,(Tf  -  Ts)/t.  (3) 

We  will  call  the  surface  thermal  forcing  given  in  the 


form  of  (3)  the  combined  boundary  condition  hereafter 
in  the  present  paper.  It  has  become  increasingly  widely 
used  as  more  reliable  data  for  Q*  has  become  available 
(Ezer  and  Mellor  1992;  Masumoto  and  Yamagata  1996; 
Mikolajewicz  and  Maier-Reimer  1994;  Sterl  and  Kat- 
tenberg  1994;  Barnier  et  al.  1995;  Stammer  et  al.  1996; 
Maltrud  et  al.  1998).  They  also  suggested  that  the  value 
of  t  can  be  interpreted  as 


pCl,Azl 
( BQ/dTs)T . 


(4) 


(Masumoto  and  Yamagata  1996;  Barnier  et  al.  1995; 
Maltrud  et  al.  1998). 

One  obvious  advantage  of  the  combined  boundary 
condition  over  the  restoring  boundary  condition  is  that, 
if  the  model  can  predict  SST  accurately,  the  heat  flux 
in  the  model  approaches  the  climatological  heat  flux. 
As  the  restoring  term  becomes  negligible,  the  model 
SST  becomes  controlled  by  the  ocean  dynamics  rather 
than  the  restoring  term  (Barnier  1998).  This  will  help 
to  remove  most  problems  associated  with  the  restoring 
boundary  condition  mentioned  above. 

This  implies,  however,  that  the  restoring  term  in  (3) 
should  be  sufficiently  small  compared  to  Q*  for  the 
combined  boundary  condition  to  really  yield  improved 
results  of  the  OGCM.  Nonetheless,  it  has  not  been  ex¬ 
amined  whether  the  restoring  term  is  sufficiently  smaller 
than  Q*  in  the  aforementioned  applications  of  the  com¬ 
bined  boundary  condition.  Of  course,  this  cannot  be 
obtained  simply  by  increasing  t  because  the  simulated 
SST  Ts  deviates  farther  from  the  climatological  SST 
7’*  as  the  restoring  becomes  weaker  with  increasing  t. 

The  reduction  of  the  restoring  term  (7’*  —  Ts)/~  can 
be  obtained  only  by  predicting  the  SST  accurately  from 
the  model.  The  most  important  factor  in  the  accurate 
prediction  of  SST  is  the  ocean  mixed  layer  process, 
which  controls  the  vertical  transfer  of  heat  below  the 
sea  surface.  This  suggests  that  the  improved  OGCM 
with  realistic  surface  heat  flux  and  sufficiently  small 
restoring  term  cannot  be  obtained  simply  by  switching 
from  the  restoring  boundary  condition  to  the  combined 
boundary  condition,  but  should  be  accompanied  by  the 
improvement  of  the  mixed  layer  process  in  the  OGCM 
as  well. 

Recently,  Stammer  et  al.  (1996),  Maltrud  et  al.  (1998), 
and  Jiang  et  al.  (1999)  compared  the  results  of  OGCMs 
using  the  restoring  and  combined  boundary  conditions 
and  observed  some  improvements  such  as  enhancement 
of  western  boundary  currents  and  the  poleward  heat 
transport.  However,  it  is  not  yet  clear  whether  the  com¬ 
bined  boundary  condition  can  provide  more  realistic 
upper-ocean  conditions,  such  as  in  the  predictions  of 
SST  and  surface  heat  flux,  than  the  restoring  boundary 
condition.  Moreover,  they  did  not  account  for  the  effects 
of  mixed  layer  processes. 

Various  ocean  mixed  layer  models  (OMLM)  have 
been  developed  to  predict  SST  and  mixed  layer  depth 
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(MLD),  usually  for  a  horizontally  homogeneous  upper 
ocean.  Since  Pacanowski  and  Philander’s  (1981)  simple 
prescription.  For  vertical  mixing  affected  by  stratifica¬ 
tion,  several  OMLMs  have  been  embedded  into  the 
OGCM,  and  their  effects  were  tested.  The  embedded 
OMLM  is  either  the  turbulence  closure  model  (Mellor 
and  Yamada  1982;  Large  et  al.  1994;  Gaspar  et  al.  1990) 
or  the  bulk  model  (Kraus  and  Turner  1967;  Price  et  al. 
1986).  The  former  is  usually  applied  to  the  z-  or  cr- 
coordinate  models,  and  the  latter  to  the  isopycnal  mod¬ 
els.  For  example,  Rosati  and  Miyakoda  (1988)  and  Ezer 
(2000)  examined  the  Mellor-Yamada  model  in  the 
OGCM,  Large  et  al.  (1997)  and  Large  and  Gent  (1999) 
examined  the  K-profile  parameterization  (KPP)  model 
by  Large  et  al.  (1994),  and  Blanke  and  Delecluse  (1993) 
examined  the  model  by  Gaspar  et  al.  (1990).  On  the 
other  hand,  Sterl  and  Kattenberg  (1994),  Oberhuber 
(1993),  and  Bleck  et  al.  (1989)  combined  bulk-type 
models,  similar  to  Kraus  and  Turner  (1967),  into  an 
OGCM.  Meanwhile,  Chen  et  al.  (1994)  suggested  an¬ 
other  bulk  model  that  hybridizes  the  Kraus-Turner  mod¬ 
el  and  the  model  by  Price  et  al.  (1986). 

Nevertheless,  it  has  not  been  shown  clearly  that  the 
prediction  of  the  global  SST  could  be  significantly  im¬ 
proved  by  including  the  OMLM.  For  example,  the  Mel¬ 
lor-Yamada  model  was  found  to  cause  excessively  high 
SST  under  heating,  and  thus  it  could  not  provide  much 
improvement  compared  to  the  case  without  the  OMLM 
under  monthly  mean  surface  forcing  (Rosati  and  Mi¬ 
yakoda  1988).  Large  et  al.  (1997)  reported  a  less  than 
5%  decrease  in  the  rms  of  the  potential  temperature 
anomaly  over  the  upper  903  m  when  the  KPP  model 
was  included,  but  they  did  not  examine  the  effects  of 
the  OMLM  on  spatial  and  temporal  variations  of  the 
SST  anomaly.  Meanwhile,  in  most  other  cases  the  ap¬ 
plication  was  restricted  to  a  regional  ocean  (Blanke  and 
Delecluse  1993;  Sterl  and  Kattenberg  1994;  Chen  et  al. 
1994;  Large  and  Gent  1999;  Ezer  2000). 

Therefore,  in  this  paper  we  attempted  to  reproduce 
more  realistic  upper  ocean  processes  such  as  SST  and 
MLD  over  the  global  ocean  by  including  a  new  OMLM 
in  an  OGCM  and  to  realize  more  realistic  surface  heat 
flux  in  the  OGCM.  In  this  way  we  hope  to  provide  an 
OGCM  that  can  be  naturally  coupled  to  an  AGCM. 

Investigation  was  made  through  three  stages  of  ex¬ 
periments.  First,  we  examined  an  OGCM  under  both 
restoring  and  combined  boundary  conditions  for  tem¬ 
perature  for  the  case  of  constant  vertical  mixing.  Sec¬ 
ond,  we  investigated  OGCMs  under  the  combined 
boundary  condition  with  different  mixed  layer  process¬ 
es:  the  constant  vertical  mixing,  the  widely  used  vertical 
mixing  scheme  by  Pacanowski  and  Philander  (1981), 
and  a  new  OMLM.  Here  a  new  OMLM  was  developed 
for  the  OGCM  based  on  the  model  by  Noh  and  Kim 
(1999).  Finally  we  repeated  the  experiments  with  longer 
restoring  timescale  for  the  case  of  the  OGCM  with  the 
new  OMLM.  The  results  were  analyzed  from  the  per¬ 
spective  of  the  improvement  of  the  upper-ocean  process. 


2.  Model  and  experiment 

a.  OGCM 


This  study  employed  the  Pacanowski  et  al.  (1991)  ver¬ 
sion  of  the  Bryan-Cox-Semtner  OGCM.  The  model  do¬ 
main  covered  the  global  ocean  up  to  80°  latitude  in  both 
hemispheres.  The  horizontal  grid  spacing  was  1°  in  both 
latitude  and  longitude.  The  vertical  grids  were  made  of  21 
levels  whose  thickness  increased  with  depth  from  20  m 
at  the  surface  to  1000  m  at  the  bottom.  To  avoid  instability 
at  high  latitudes  a  simple  symmetric  finite  impulse  filter 
was  applied  in  both  hemispheres  beyond  55°.  Realistic 
bottom  topography  was  used  as  much  as  possible  within 
the  limits  of  grid  resolution,  including  nine  islands  and  a 
rather  realistic  Indonesian  archipelago. 

Horizontal  eddy  viscosity  AM  and  diffusivity  AH  were 
1  X  10s  cm2  s_I  and  1  X  107  cm2  s-1,  respectively. 
Fluxes  of  heat  and  salinity  were  fixed  to  zero  at  the 
bottom  and  at  the  side  wall.  Bottom  friction  was  given 
by  CD  lulu  with  CD  =  0.0025,  and  a  no-slip  boundary 
condition  was  applied  at  the  sidewall. 

Vertical  eddy  viscosity  KM  and  diffusivity  KH  were 
given  by  1  and  0.2  cm2  s'1,  respectively,  in  the  case  of 
constant  vertical  mixing.  Convection  is  implemented  by 
applying  very  large  eddy  diffusivity  (KH  =  105  cm2 
s  1 ),  whenever  unstable. 

In  the  case  with  the  Pacanowski  and  Philander  scheme 
(1981)  Km  and  K  H  were  given  by 


Km  (1  +  oRi)2  +  Kmo ’ 

K » -  rns + 


(5) 

(6) 


where  the  empirical  constants  were  given  by  a  =  5,  K0 
=  50  cm2  s_1,  Km0  =  0.0134  cm2  s_1,  and  Km  =  0.001 
34  cm2  s  1 .  Here  Ri  is  the  Richardson  number 

N2 

Ri  =  - ,  (7) 

(■ du/dz Y  +  (dvldz)2 

where  N  is  the  Brunt-Vaisala  frequency. 

The  model  was  started  from  a  state  of  rest,  with  an¬ 
nual  mean  temperature  and  salinity  distributions  by  Lev- 
itus  (1982),  and  forced  by  Hellerman  and  Rosenstein 
(1983)  wind  stress  climatology.  To  prescribe  the  heat 
flux  at  the  surface  the  reanalysis  data  from  the  National 
Centers  for  Environmental  Prediction  (NCEP)  were 
used.  Meanwhile,  optimally  interpolated  SST  (OISST) 
(Reynolds  and  Smith  1994)  data  were  used  for  SST 
instead  of  Levitus  data  for  more  accurate  prescription 
of  the  surface  boundary  condition.  For  the  surface  forc¬ 
ing  the  monthly  mean  data  were  used  with  a  linear  in¬ 
terpolation  at  each  time  step.  Heat  flux,  given  by  Q  in 
(3),  was  set  to  zero  when  the  SST  is  below  — 2°C,  as¬ 
suming  that  the  sea  surface  is  covered  by  ice.  Surface 
salinity  was  restored  to  the  climatological  value  of  Lev¬ 
itus  (1982)  with  a  restoring  timescale  of  30  days 
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(a)  EXP  AO 


60°E  120°E  180"  120°W  60°W 


(b)  EXP  A1 
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Fig.  1.  Annual  mean  volume  transport  streamfunction  (Sv):  (a)  EXP  AO  and  (b)  EXP  Al. 


throughout  all  experiments.  The  model  was  integrated 
for  24  years,  which  is  believed  to  be  sufficient  for  ex¬ 
amining  the  changes  in  the  heat  flux  and  temperature 
in  the  upper  ocean  (Sterl  and  Kattenberg  1994). 

b.  OMLM 

Recent  measurements  of  the  microstructure  of  the  oce¬ 
anic  boundary  layer  from  various  field  experiments  have 
revealed  that  dynamic  processes  of  the  oceanic  boundary 
layer  is  fundamentally  different  from  that  of  the  atmo¬ 
spheric  boundary  layer  (see,  e.g.,  Drennan  et  al.  1996; 
Terray  et  al.  1996).  Particularly  they  observed  that  tur¬ 
bulent  kinetic  energy  (TKE)  is  enhanced  remarkably  near 
the  sea  surface,  unlike  near  the  rigid  boundary  of  the 
atmospheric  boundary  layer.  A  high  level  of  TKE  near  the 
sea  surface  helps  to  maintain  the  well-mixed  layer  in  the 
upper  ocean.  This  also  makes  downward  flux  of  TKE 
important  in  the  TKE  budget  of  the  upper  mixed  layer. 


while  rendering  TKE  production  by  the  mean  velocity 
shear  relatively  unimportant  (Noh  1996). 

Noh  and  Kim  (1999)  recently  developed  a  new  OMLM, 
taking  these  observational  evidences  into  consideration. 
The  model  is  a  second-order  turbulence  closure  model 
using  eddy  diffusivity,  but  it  produces  a  well-mixed  layer 
even  under  the  stabilizing  heat  flux,  consistent  with  the 
observation  and  the  bulk-type  OMLMs  (i.e.,  Kraus  and 
Turner  1977),  whereas  most  other  OMLMs  using  eddy 
diffusivity  lead  to  strong  stratification  and  shear  near  the 
surface  (i.e.,  Mellor  and  Yamada  1982). 

The  model  calculates  the  vertical  eddy  viscosity  KM, 
the  eddy  diffusivity  of  temperature  and  salinity  KH,  and 
the  eddy  diffusivity  of  turbulent  kinetic  energy  KE  by 


KM  = 

sMql 

(8) 

KH  = 

SHql, 

(9) 

KE  = 

SEql, 

(10) 
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Fig.  2.  Distributions  of  the  SST  anomaly  (°C)  for  EXP  AO:  (a)  annual  mean,  (b)  Jan,  and  (c) 
Jul.  Contour  interval  is  1°C,  and  solid  and  dashed  lines  represent  positive  and  negative  anomalies, 
respectively.  Areas  with  an  SST  anomaly  larger  than  1°C  are  shaded. 
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Fig.  3.  As  in  Fig.  2  but  for  EXP  Al. 
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Fig.  4.  Seasonal  variations  of  the  mean  SST  (°C)  (climatology, 
EXP  AO,  EXP  Al):  (a)  Northern  Hemisphere,  (b)  Southern  Hemi¬ 
sphere,  and  (c)  equatorial  region  (5°S-5°N).  OISST  represents  Op¬ 
timally  Interpolated  SST  data  (Reynolds  and  Smith  1994). 


where  q2l2(=E)  is  TKE,  l  is  the  length  scale  of  tur¬ 
bulence,  and  SM,  SH,  and  SE  are  empirical  constants. 
Turbulent  kinetic  energy  E  is  obtained  from 
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Fig.  5.  As  in  Fig.  4  but  for  mean  heat  flux  (W  m  2). 


and  Yamada  model  in  the  parameterizations  of  the  sur¬ 
face  forcing  and  of  the  effects  of  stratification.  The 
boundary  condition  for  TKE  is  given  by 

dE 

Ke—  =  mu%,  (13) 

dz 


where  u *  is  the  frictional  velocity  due  to  wind  stress. 
Meanwhile,  the  length  scale  /  is  given  by 


k(z  +  Zp) 

1  +  k(z  +  z,0)/h' 


(14) 


if  horizontal  advection  and  diffusion  are  neglected.  Here 
g,  p0,  and  p  represent  the  gravitational  acceleration  and 
the  reference  and  local  density,  respectively.  The  dis¬ 
sipation  rate  e  is  calculated  by 

e  =  CBqVl.  (12) 

Equations  (8)— (12)  are  typical  for  the  second-order  tur¬ 
bulence  closure  model.  The  coefficients  are  taken  as  the 
same  as  in  Mellor  and  Yamada  (1982)  in  a  homogeneous 
fluid:  SM  =  0.39(=SMO),  SmISh  =  0.8,  and  SJSE  =  1.95, 
and  Ce  =  0.06(=CeO). 

However,  the  model  is  quite  different  from  the  Mellor 


where  h  is  the  depth  of  a  mixed  layer,  k  the  von  Karman 
constant  (  =  0.4),  and  z0  the  roughness  length  scale  at 
the  surface. 

Based  on  Craig  and  Banner’s  analysis  (1994)  on  re¬ 
cent  observational  data  of  the  upper  ocean,  Noh  and 
Kim  (1999)  assumed  the  values  m  =  100  and  z0  =  1 
m  for  (13)  and  (14).  These  values  characterize  the  in¬ 
tensive  mixing  near  the  sea  surface,  as  mentioned  above. 
In  this  case  the  relevant  Richardson  number,  which  is 
used  to  represent  the  effects  of  stratification,  should  be 
determined  in  terms  of  TKE  itself  rather  than  the  mean 
velocity  shear  because  the  TKE  flux  plays  an  important 
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Fig.  6.  Meridional  overturning  streamfunction  (Sv):  (a)  EXP  AO 
and  (b)  EXP  Al. 

role  in  TKE  production.  Therefore,  the  TKE  Richardson 
number  defined  by 

R,  =  c Nl/qy-  (15) 

must  be  used  instead  of  Ri  given  by  (7). 

In  this  case  we  can  parameterize  the  effects  of  strat¬ 
ification  as 


SM/SM0  =  (1  +  «R,)-"2, 

(16) 

CJCB0  =  (1  +  aR,)-U2, 

(17) 

similar  to  Wyngaard  (1985),  where  a  is  an  empirical 
constant.  Meanwhile,  the  proportionality  of  SH  and  .S', 
to  SM  is  maintained. 

A  uniform  mixed  layer  is  maintained  during  convec¬ 
tion,  which  means  KM  =  KH  =  °°  whenever  R,  <  0. 
However,  the  vertical  transfer  of  TKE  is  assumed  to  be 
unaffected  by  convection,  that  is,  KE  =  KE(R,  =  0), 
based  on  the  scaling  analysis.  See  Noh  and  Kim  (1999) 
for  more  detailed  description  of  the  model. 

The  resultant  model  was  successful  for  the  first  time 
in  predicting  correctly  the  evolution  of  both  the  vertical 
profiles  of  temperature  and  the  turbulence  structure  of 
the  mixed  layer,  as  well  as  the  SST,  that  was  the  main 
concern  in  the  previous  results  (Noh  and  Kim  1999). 


The  new  OMLM  was  also  shown  to  realistically  sim¬ 
ulate  the  formation  of  a  diurnal  thermocline,  which  re¬ 
mains  one  of  the  most  difficult  problems  in  most 
OMLMs  (Noh  1996).  However,  the  model  was  devel¬ 
oped  mainly  for  prediction  of  the  diurnal  variation  of 
the  upper  ocean.  Therefore,  some  modifications  are  nec¬ 
essary  so  as  to  be  applied  to  the  seasonal  variation  of 
the  deeper  ocean  as  in  the  OGCM. 

First,  we  expect  that  vertical  mixing  during  convec¬ 
tion  into  the  deeper  ocean  cannot  occur  instantaneously 
but  over  a  certain  timescale  with  the  corresponding  finite 
values  of  KH.  We  assume  the  eddy  diffusivity  during 
convection  as  KH  =  105  cm2  s-1,  following  the  sug¬ 
gestion  by  Klinger  et  al.  (1996).  It  was  reported,  how¬ 
ever,  that  the  results  are  insensitive  to  the  value  of  KH 
during  convection  as  long  as  they  are  within  the  range 
of  104-107  cm2  s_1  (Hirst  and  Cai  1994). 

Meanwhile,  Killworth  (1989)  suggested  that  it  is  not 
desirable  to  force  the  vertical  mixing  of  momentum  dur¬ 
ing  convection  in  the  OGCM  since  geostrophic  adjust¬ 
ment  tends  to  make  the  velocity  field  vertically  uniform 
for  large-scale  flows.  So  we  assume  that  eddy  viscosity 
Km,  as  well  as  KE,  maintains  the  value  of  a  homogeneous 
fluid,  unaffected  by  convection.  It  was  also  observed, 
however,  by  Killworth  that  the  results  are  largely  in¬ 
dependent  of  whether  momentum  is  mixed  or  not  during 
convection. 

Noh  and  Kim  (1999)  evaluated  the  empirical  constant 
a  as  a  —  120  in  (16)  and  (17).  However,  the  vertical 
resolution  in  that  case  was  1  m,  whereas  it  is  larger  than 
20  m  in  the  present  case.  We  expect  that  the  determi¬ 
nation  of  the  values  of  N,  and  thus  R ,  (and  also  Ri),  is 
affected  by  the  vertical  resolution  of  the  model.  This 
will  be  particularly  true  near  the  thermocline  where  the 
effects  of  stratification  play  a  critical  role.  This  makes 
it  necessary  to  reevaluate  a.  The  appropriate  value  of 
a  was  estimated  to  be  a  —  1000  from  the  analysis  of 
the  present  OGCM  results. 

Turbulence  generated  at  the  surface  cannot  be  trans¬ 
ferred  below  the  mixed  layer  (Noh  1996),  and  the  sourc¬ 
es  of  turbulence  below  the  mixed  layer  are  fundamen¬ 
tally  different  from  those  within  the  mixed  layer.  There 
are  various  sources  of  mixing  below  the  mixed  layer: 
for  example,  internal  wave  breaking,  vertical  shear,  and 
double  diffusion.  Although  various  parameterizations 
could  be  suggested  to  represent  the  mixing  owing  to  the 
afore-mentioned  mechanisms  (Large  et  al.  1994;  Toole 
1998),  there  still  remains  extreme  uncertainty  regarding 
how  important  they  are  and  how  they  should  be  param¬ 
eterized.  Moreover,  the  previously  mentioned  sources 
of  turbulence  are  related  to  small-scale  phenomena  that 
require  fine-scale  resolution  for  meaningful  parameter¬ 
ization,  whereas  the  vertical  resolution  of  the  OGCM  is 
usually  very  poor  in  the  deep  ocean,  having  a  typical 
vertical  grid  size  of  a  few  hundred  meters. 

Therefore,  in  this  model  we  assume  that  the  eddy 
viscosity  KMB  and  eddy  diffusivity  KHB  below  the  mixed 
layer  can  be  roughly  estimated  as  1  and  0.1  cm2  s  1 , 
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Fig.  7.  Distributions  of  the  SST  anomaly  (°C)  for  EXP  Bl:  (a)  annual  mean,  (b)  Jan,  and  (c) 
Jul.  The  contour  interval  is  1°C,  and  solid  and  dashed  lines  represent  positive  and  negative 
anomalies,  respectively.  Areas  with  an  SST  anomaly  larger  than  1°C  are  shaded. 
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Fig.  8.  As  in  Fig.  7  but  for  EXP  Cl. 
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Fig.  9.  Seasonal  variations  of  the  mean  SST  (°C)  (climatology, 
EXP  Al,  EXP  Bl,  EXP  Cl):  (a)  Northern  Hemisphere,  (b)  Southern 
Hemisphere,  and  (c)  equatorial  region  (5°S-5°N). 


Fig.  10.  As  in  Fig.  9  but  for  the  rms  of  SST  anomaly  (°C)  (EXP 
Al,  EXP  Bl,  EXP  Cl). 


respectively,  as  the  simplest  way  of  representing  sub¬ 
surface  mixing  under  this  extremely  uncertain  situation. 
These  values  follow  the  estimation  of  the  mixing  due 
to  internal  wave  breaking  given  by  Large  et  al.  (1994) 
and  Toole  (1998),  estimated  from  the  observational  data 
(Gregg  1989;  Peters  et  al.  1988;  Ledwell  et  al.  1993) 
and  the  Garrett-Munk  (1975)  internal  wave  model.  On 
the  other  hand,  the  magnitudes  of  mixing  owing  to  ver¬ 
tical  shear  and  double  diffusion  are  expected  to  be  not 
larger  than  that  by  internal  wave  breaking,  although  it 
is  extremely  difficult  to  estimate  them  (Toole  1998).  It 
is  thus  assumed  that  the  contributions  to  mixing  from 
these  mechanisms  are  included  in  the  above  estimations 
of  Kmb  and  KHB. 

Noh  (1996)  suggested  that  the  mixed  layer  depth  can 
be  represented  most  appropriately  by  the  depth  where 
the  vertical  mixing  decreases  to  a  certain  level.  Ac¬ 
cordingly,  in  the  present  model  the  definition  of  mixed 
layer  depth  h  is  given  by  the  depth  at  which  the  cal¬ 
culated  eddy  viscosity  KH  decreases  with  depth  to  the 
value  of  Khb.  For  depths  deeper  than  li,  we  simply  put 
Km  —  Kmb  and  K H  —  K HB. 


As  described  above,  the  present  OMLM  actually  pre¬ 
scribes  the  vertical  mixing  over  the  whole  depth  of  the 
ocean  as  well  as  within  the  mixed  layer,  although  most 
concerns  focus  on  the  latter.  Therefore,  the  present 
OMLM  can  be  also  referred  to  as  the  vertical  mixing 
scheme,  as  are  the  cases  of  other  OMLMs  (Rosati  and 
Miyakoda  1988;  Large  et  al.  1997). 


c.  Experiments 

Numerical  experiments  with  different  combinations 
of  the  mixed  layer  process  and  surface  thermal  forcing 
were  carried  out  in  the  following  three  steps.  First,  we 
considered  the  effects  of  different  surface  boundary  con- 


Table  1.  Rms  difference  of  the  simulated  regional  mean  SST  (C°) 
from  the  observed  data  in  the  Northern  Hemisphere  (NH),  the  South¬ 
ern  Hemisphere  (SH),  and  the  equatorial  region  (EQ). 


Expt 

NH 

SH 

EQ 

EXP  Al 

0.45 

0.69 

0.83 

EXP  Bl 

0.32 

0.46 

0.64 

EXP  Cl 

0.24 

0.17 

0.34 
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Table  2.  As  in  Table  1  but  for  the  regional  mean  of  the  rms 
difference  of  the  simulated  SST  (C°)  from  the  observed  data. 


Expt 

NH 

SH 

EQ 

EXP  Al 

1.37 

1.29 

1.90 

EXP  B1 

1.39 

1.22 

1.88 

EXP  Cl 

1.28 

1.03 

1.50 

ditions  for  temperature  using  the  OGCM  with  constant 
vertical  mixing  coefficients,  that  is,  the  restoring  bound¬ 
ary  condition  (EXP  AO)  and  the  combined  boundary 
condition  (EXP  Al).  In  both  cases  the  restoring  time- 
scale  rwas  fixed  to  30  days.  The  results  were  analyzed 
from  the  perspective  of  the  realizations  of  the  observed 
SST  and  heat  flux. 

Second,  we  carried  out  the  experiments  with  OGCMs 
under  the  same  combined  boundary  condition  but  with 
different  vertical  mixing  schemes.  We  tried  the  vertical 
mixing  schemes  by  Pacanowski  and  Philander  (1981) 
and  the  new  OMLM  explained  in  the  previous  section, 
referred  to  as  EXP  B1  and  EXP  Cl,  respectively.  By 
comparing  the  results  from  EXP  B1  and  EXP  Cl  to¬ 
gether  with  those  from  EXP  Al  we  investigated  the 
effects  of  the  mixed  layer  process  on  the  results  of  the 
OGCM,  with  major  emphasis  on  the  upper  ocean. 

Note  that,  when  the  restoring  boundary  condition  is 
used  in  the  OGCM,  the  more  accurate  prediction  of  SST 
means  the  poorer  prediction  of  surface  heat  flux,  which 
makes  it  impossible  to  compare  the  performance  of  the 
OMLM.  On  the  other  hand,  in  the  combined  boundary 
condition  the  prediction  of  more  accurate  SST  coincides 
with  more  realistic  surface  heat  flux,  and  thus  illustrates 
the  better  performance  of  the  OGCM. 

Finally,  once  we  realized  that  the  new  OMLM  per¬ 
forms  well,  we  repeated  EXP  Cl  but  with  a  restoring 
timescale  increased  to  a  more  realistic  value  (t  =  60 
days),  referred  to  as  EXP  C2.  From  the  comparison 
between  EXP  Cl  and  EXP  C2  we  investigated  how  the 
increased  restoring  time  affects  the  general  performance 
of  the  model.  Note  that  the  OMLM  should  perform  well 
under  the  natural  feedback  between  the  atmosphere  and 
the  ocean  if  it  really  represents  the  mixed  layer  process 
properly. 

3.  Results 

a.  Comparison  between  the  restoring  and  combined 

surface  thermal  forcing  (EXPs  AO  and  Al ) 

Figure  1  shows  the  annual  mean  volume  transport 
streamfunction  obtained  from  EXPs  AO  and  Al,  re¬ 
spectively.  Both  cases  are  similar  to  each  other  and  re¬ 
produce  the  general  features  of  the  global  circulation  of 
the  ocean  satisfactorily  regardless  of  the  difference  in 
the  surface  boundary  condition  of  temperature.  This 
confirms  that  the  global  circulation  is  driven  primarily 
by  wind  stress.  In  both  cases  the  intensity  of  western 
boundary  currents  are  somewhat  weaker  and  their  sep- 
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Fig.  11.  Meridional  overturning  streamfunction  (Sv):  (a)  EXP  B1 
and  (b)  EXP  Cl. 

aration  points  are  located  rather  northward  compared  to 
the  observation,  probably  due  to  insufficient  grid  res¬ 
olution,  as  in  the  cases  of  most  other  OGCMs  (see,  e.g., 
Haidvogel  and  Bryan  1992).  Nevertheless,  the  maxi¬ 
mum  annual  mean  transports  of  the  Gulf  Stream  and 
the  Kuroshio,  given  by  approximately  30-40  Sv  and 
50-60  Sv  (Sv  =  106  m3  s~')  are  in  agreement  with  the 
high-resolution  simulation  results  (Semtner  and  Chervin 
1992)  and  other  observation  data  (Tomczak  and  Godfrey 
1994).  The  transports  of  the  Antarctic  Circumpolar  Cur¬ 
rent  (ACC)  through  Drake  Passage  and  the  Indonesian 
Throughflow  are  approximately  100-110  Sv  and  10- 
20  Sv,  respectively,  and  are  also  within  the  reasonable 
range  of  observational  estimates  (Fieux  et  al.  1994; 
Godfrey  et  al.  1993;  Nowlin  and  Klinck  1986). 

The  only  noticeable  difference  is  observed  in  the 
high-latitude  region  where  deep  convection  occurs.  The 
convective  pattern  and  the  consequential  anticlockwise 
circulation  in  the  northern  Atlantic  appear  in  EXP  Al, 
but  not  in  EXP  A0,  which  may  be  associated  with  the 
fact  that  convection  is  suppressed  by  the  restoring 
boundary  condition  (Mikolajewicz  and  Maier-Reimer 
1994;  Cai  and  Chu  1996;  Toggweiler  et  al.  1989).  It  is 
also  observed  that  the  ACC  is  slightly  stronger  in  EXP 
A0,  whereas  the  Gulf  Stream  is  slightly  stronger  in  EXP 
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Fig.  12.  Cross  sections  of  the  annual-mean  zonal  velocity  at  the  equator  and  155°W:  (a) 
climatology  [Hawaii-Tahiti  Shuttle  Experiment  (Wyrtki  and  Kilonsky  1984)],  (b)  EXP  Al,  (c) 
EXP  Bl,  and  (d)  EXP  Cl. 


Al.  The  differences  in  OGCM  results  between  two  dif¬ 
ferent  surface  forcing  methods  mentioned  above  are  also 
in  agreement  with  the  results  by  Maltrud  et  al.  (1998), 
who  used  higher-resolution  grids  and  higher-frequency 
data. 

It  should  be  mentioned,  however,  that  the  circulation 
of  the  deep  ocean  is  still  controlled  by  the  initial  T,  S 
distributions  given  by  the  observed  data  in  the  present 
simulation.  Therefore,  we  expect  that  the  volume  trans¬ 
port  streamfunction  of  Fig.  1  and  the  meridional  over¬ 
turning  streamfunction,  given  below  by  Fig.  5,  are  still 
affected  by  the  initial  T,  S  distributions.  Nevertheless, 
they  still  serve  effectively  to  demonstrate  the  effects  of 
surface  forcing  and  the  mixed  layer  process  because  the 
circulations  of  the  deep  ocean  due  to  the  initial  T,  S 
distributions  remain  equivalent  in  any  case. 

Figures  2  and  3  show  the  global  distributions  of  SST 
anomaly  from  the  observed  distribution,  for  the  cases 
of  the  annual  mean,  January  and  July,  from  EXPs  AO 
and  Al.  The  significant  underestimation  of  SST  in  the 
equatorial  ocean  is  observed  in  both  cases.  This  is  very 
common  in  most  OGCMs  and  is  due  to  the  too  strong 
South  Equatorial  Current  (SEC)  and  upwelling.  We  will 
discuss  this  again  in  the  next  section.  The  underesti¬ 
mation  of  SST  in  the  equatorial  ocean  is  somewhat  im¬ 
proved  in  EXP  Al,  mainly  because  larger  heat  flux  is 
applied  in  that  case  (see  also  Fig.  4). 

Significant  SST  anomalies  are  also  found  in  both  cas¬ 
es  in  the  regions  where  the  simulated  circulation  de¬ 
viates  from  the  real  circulation,  such  as  the  separation 
points  of  the  western  boundary  currents.  For  example, 
the  simulated  Gulf  Stream,  which  is  too  broad  and  over¬ 
shoots  beyond  the  separation  point,  causes  the  dipole 
SST  anomaly,  as  pointed  out  by  Sterl  and  Kattenberg 
(1994).  That  is,  water  along  the  coast  appears  too  warm 
and  recirculated  water  over  the  region  of  warm  water 


moving  northward  too  cool.  The  displacement  of  sep¬ 
aration  points  of  the  Kuroshio,  Brazil,  and  Agulhas  Cur¬ 
rents  also  causes  the  large  SST  anomalies  in  the  regions 
off  Japan,  Argentina,  and  South  Africa.  Meanwhile,  a 
too  weak  Gulf  Stream,  which  does  not  extend  to  the 
North  Sea,  causes  the  negative  SST  anomaly  in  the 
North  Sea.  A  large  SST  anomaly  in  the  ACC  region 
implies  an  inaccurate  representation  of  variability  of  the 
ACC  and  Polar  Front.  Hopefully,  these  shortcomings 
can  be  remedied  by  more  accurate  simulation  of  cir¬ 
culation,  possibly  with  a  higher-resolution  OGCM. 

However,  the  most  striking  SST  anomaly  is  observed 
in  the  seasonal  variation  of  SST,  where  the  opposite 
tendency  of  the  SST  anomaly  appears  between  two  cas¬ 
es.  The  SST  turns  out  to  be  too  low  during  summer  and 
too  high  during  winter  in  EXP  AO  (Figs.  2b,c),  and  vice 
versa  in  EXP  Al  (Figs.  3b, c).  This  means  that  the  am¬ 
plitude  of  the  seasonal  cycle  of  SST  is  too  small  in  EXP 
AO,  whereas  it  is  too  large  in  EXP  Al,  as  shown  in  the 
seasonal  variations  of  the  hemispheric  mean  SST  (Fig. 

4) . 

In  the  case  of  the  restoring  boundary  condition  (EXP 
AO),  the  SST  is  actually  “dragged”  to  the  climatological 
SST  and  thus  results  in  the  reduced  amplitude  and  phase 
lag  of  the  seasonal  variation  of  SST,  as  mentioned  by 
Haidvogel  and  Bryan  (1993).  On  the  other  hand,  over¬ 
estimation  of  SST  during  the  warming  season  in  EXP  Al 
is  caused  by  improper  vertical  transfer  of  heat  from  the 
surface.  Although  we  can  attempt  to  improve  the  predic¬ 
tion  of  SST  by  increasing  the  vertical  mixing  coefficient 
in  this  case,  this  leads  to  a  severely  unrealistic  vertical 
structure  of  temperature,  which  already  suffers  from  a  too 
deep  and  diffuse  thermocline  (Haidvogel  and  Bryan  1993). 
During  summer  too  high  SST  in  EXP  Al  causes  the  heat 
flux  to  be  underestimated,  as  expected  from  (3)  (see  Fig. 

5) .  Consequently,  an  insufficient  amount  of  heat  accu- 
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latitude 

Fig.  12.  (Continued) 


mulated  during  summer  possibly  leads  to  too  low  SST 
under  the  cooling  in  winter  in  Fig.  3. 

Figure  5  shows  the  seasonal  cycle  of  the  resultant 
hemispheric  mean  heat  flux  based  on  the  simulated  SST 
in  comparison  with  the  climatological  heat  flux  Q*.  The 
resultant  heat  flux  is  found  to  be  improved  in  EXP  Al, 
as  expected  from  the  fact  that  the  climatological  heat 
flux  itself  is  included  in  the  calculation.  Nonetheless, 
the  contribution  from  the  restoring  term  still  remains 
significantly  large,  and  this  suggests  that  many  problems 
associated  with  the  restoring  boundary  condition  still 
persist.  The  experiment  with  the  larger  restoring  time- 
scale  of  60  days  aggravates  the  situation  further  (not 


shown)  because  the  simulated  SST  deviates  severely 
from  the  observed  SST. 

As  expected  from  Fig.  5,  the  increased  seasonal  cycles 
of  heat  flux  result  in  a  stronger  meridional  overturning 
streamfunction  (Fig.  6).  Especially  the  suppression  of 
deep  convection  in  the  North  Atlantic,  as  mentioned 
above,  is  evidenced  clearly  in  Fig.  6a  and  it  causes  the 
clockwise  circulation  in  the  northern  thermohaline  cell 
at  intermediate  depths  to  be  very  weak  in  EXP  AO.  The 
situation  is  improved  in  EXP  Al.  Meanwhile,  the  over¬ 
all  results  in  Fig.  6  are  in  good  agreement  with  Semtner 
and  Chervin  (1992)  and  Maltrud  et  al.  (1998),  but  the 
thermohaline  circulation  in  the  Northern  Hemisphere  is 
weaker  than  that  in  Large  et  al.  (1997),  where  an  ac¬ 
celerated  integration  was  applied  to  achieve  equilibrium 
in  the  deep  ocean. 

b.  Effects  of  the  mixed  layer  process  (EXPs  Al,  Bl, 

and  Cl ) 

In  this  section  we  compare  the  effects  of  the  mixed 
layer  process  based  on  the  results  from  OGCMs  with 
constant  vertical  eddy  diffusivity  and  viscosity  (EXP 
Al),  with  the  Pacanowski  and  Philander  (PP)  scheme 
(EXP  B 1),  and  with  the  new  OMLM  described  in  section 
2  (EXP  Cl). 

Figures  7  and  8  show  the  global  distributions  of  SST 
anomalies  from  the  observed  distribution  for  the  cases 
of  the  annual  mean,  January,  and  July  from  EXPs  Bl 
and  Cl.  They  clearly  show,  together  with  Fig.  3,  that 
the  seasonal  variation  of  SST  anomaly  is  significantly 
affected  by  the  vertical  mixing  process.  The  most  re¬ 
markable  improvement  due  to  including  the  new  OMLM 
is  the  disappearance  of  the  strong  seasonal  bias  observed 
in  EXP  Al. 

Although  considerable  improvement  is  obtained  in 
EXP  Bl  (Fig.  7),  the  results  from  the  PP  scheme  still 
show  an  overestimated  amplitude  of  SST  cycle.  In  par¬ 
ticular,  they  give  poor  results  near  the  ACC  where  the 
wind  stress  is  strong  (Trenberth  et  al.  1989),  thus  leading 
to  poor  performance  in  the  SST  cycle  in  the  Southern 
Hemisphere.  Note  that  the  PP  scheme  does  not  account 
for  the  effects  of  wind  stress  on  vertical  mixing.  Also, 
the  PP  scheme  causes  a  rather  large  SST  anomaly  in 
the  eastern  Pacific  where  SST  is  known  to  be  affected 
by  the  seasonally  varying  wind  stress  (Chen  et  al.  1994). 

On  the  other  hand,  it  is  observed  that  inclusion  of  the 
new  OMLM  increases  slightly  the  SST  anomaly  in  re¬ 
gions  where  the  simulated  circulation  deviates  from  the 
real  circulation,  such  as  the  regions  of  the  western 
boundary  current  separation.  With  the  inclusion  of 
OMLM  the  values  of  KM  and  K H  vary  in  space  and  time, 
and  therefore  the  inaccurate  prediction  of  the  ocean  cir¬ 
culation  may  exacerbate  the  situation. 

In  order  to  represent  the  performance  of  different  ver¬ 
tical  mixing  schemes  in  the  OGCM  more  quantitatively 
we  obtained  both  seasonal  variations  of  the  mean  SST 
(Fig.  9)  and  the  rms  SST  anomaly  (Fig.  10)  for  both 
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Fig.  13.  Distributions  of  the  MLD  (m)  (climatology).  MLD  is  defined  by  the  depth  of  the  temperature 
difference  1°C  from  SST:  (a)  Jan  and  (b)  Jul. 


hemispheres  and  the  equatorial  region.  The  rms  SST 
anomaly  may  represent  the  actual  performance  in  the 
present  OGCM  (Fig.  10),  but  these  values  also  include 
the  error  due  to  the  inaccurate  circulation,  as  mentioned 
in  the  previous  section.  On  the  other  hand,  the  general 
performance  of  vertical  mixing  scheme,  that  is,  how  the 
SST  over  the  global  ocean  responds  to  the  surface  heat¬ 
ing  and  cooling  in  the  one  dimensional  sense,  may  be 
judged  by  seasonal  variations  of  the  mean  SST  (Fig.  9). 
Both  Figs.  9  and  10  clearly  show  that  the  prediction  of 
SST  could  be  improved  significantly  by  introducing  a 
new  OMLM  in  all  regions. 

To  summarize  the  results  shown  in  Fig.  9,  the  rms 
difference  of  the  regional  mean  SST  from  the  obser¬ 


vation  data,  which  is  given  by  [£F ,  (Tt  —  Tf)2/ 1 2] 1/2 , 
where  the  overbar  represents  the  regional  mean  and  i 
represents  a  month,  was  calculated  and  the  results  are 
listed  in  Table  1.  Similarly  annual  mean  values  of  the 
rms  difference  of  the  simulated  SST  from  the  observed 
data,  as  shown  in  Fig.  10,  are  listed  in  Table  2.  These 
data  again  corroborate  the  improvement  of  the  perfor¬ 
mance  with  the  new  OMLM.  Note  that  the  performance 
of  the  rms  difference  of  SST  is  not  as  good  as  that  of 
the  mean  SST,  because  it  includes  the  error  owing  to 
the  inaccurate  circulation.  Moreover,  the  effect  of  the 
inaccurate  circulation  on  the  SST  anomaly  becomes 
larger  in  the  presence  of  the  OMLM,  as  mentioned  ear¬ 
lier. 
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Fig.  14.  As  in  Fig.  13  but  for  EXP  Al. 
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The  substantial  improvement  of  the  annual  mean  and 
seasonal  cycle  of  SST  with  the  new  OMLM  in  the  equa¬ 
torial  region  is  also  observed  in  Figs.  8,  9,  and  10.  When 
an  experiment  was  carried  out  with  the  wind  stress  of 
Hellerman  and  Rosenstein  (1983)  reduced  to  0.75,  as 
suggested  by  Stockdale  et  al.  (1993)  for  a  better  pre¬ 
scription  of  the  wind  stress  in  the  equatorial  region,  it 
was  observed  that  the  SST  prediction  in  the  equatorial 
region  is  improved  even  further  (not  shown). 

This  improvement  in  the  equatorial  region  is  attri¬ 
buted  to  the  realization  of  a  more  realistic  mixed  layer 
in  the  equatorial  region.  The  mixed  layer  with  appro¬ 
priate  thickness  prevents  an  accumulation  of  momentum 
and  the  consequent  overestimation  of  the  equatorial  cur¬ 
rents  by  distributing  the  momentum  vertically.  This 


helps  to  suppress  the  excessive  equatorial  upwelling.  It 
can  be  clearly  understood  from  the  comparison  of  the 
meridional  streamfunctions  (Figs.  6  and  11).  In  the  case 
of  EXP  Al  the  equatorial  upwelling  is  much  stronger, 
and  the  equatorial  divergence  is  limited  to  the  abnor¬ 
mally  shallow  surface  layer  in  the  equatorial  region. 
Whereas,  EXP  Cl  shows  weaker  upwelling  and  equa¬ 
torial  divergence  over  a  thickness  of  about  100  m,  which 
is  comparable  to  the  observed  mixed  layer  depth  (Hal- 
perin  1980).  The  reduced  upwelling  due  to  including 
the  mixed  layer  process  was  also  observed  by  Sterl  and 
Kattenberg  (1994).  Also,  it  is  observed  that  the  depth 
to  which  equatorial  upwelling  occurs  becomes  shallower 
in  EXP  Cl  than  in  EXP  Al.  The  PP  scheme  shows  the 
intermediate  pattern  in  intensity  of  equatorial  diver- 
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Fig.  15.  As  in  Fig.  13  but  for  EXP  Bl. 


gence,  but  the  vertical  extension  of  the  equatorial  up- 
welling  is  even  smaller  than  in  EXP  Cl. 

To  examine  the  effects  of  vertical  mixing  in  the  equa¬ 
torial  region  in  more  detail,  we  plotted  the  zonal  velocity 
distributions  at  about  155°W  at  the  equator  for  three 
cases  (Fig.  12).  The  observation  data  shows  that  the 
Equatorial  Undercurrent  (EUC)  is  generally  limited  be¬ 
tween  50  and  250  m  with  maximum  velocity  about  90 
cm  s_1  (Wyrtki  and  Kilonsky  1984).  In  the  case  of  EXP 
A1  the  EUC  is  vertically  exaggerated  and  extends  very 
close  to  the  sea  surface,  which  is  typical  in  the  OGCM 
without  the  mixed  layer  process  (Semtner  and  Chervin 
1992).  On  the  other  hand,  in  the  case  of  EXP  Cl,  the 
depth  of  the  EUC  is  generally  comparable  to  observa¬ 
tions,  although  it  tends  to  be  slightly  deeper  than  ob¬ 


served.  Downward  overextension  of  the  EUC  to  a  depth 
of  about  350  m  in  both  EXP  A1  and  EXP  Cl  is  due  to 
the  low  vertical  resolution  of  the  present  model  having 
a  grid  thickness  of  100  m  at  that  depth.  The  low  vertical 
resolution  may  also  cause  the  underestimation  of  the 
intensity  of  the  EUC.  The  case  of  EXP  Bl  is  similar  to 
that  of  EXP  Cl,  but  the  EUC  is  observed  to  extend 
downward  too  deep  compared  to  the  observation,  sug¬ 
gesting  improper  vertical  mixing  below  the  EUC. 

Annual-mean  volume  transport  streamfunctions  from 
both  EXP  Bl  and  EXP  Cl  are  similar  to  that  from  EXP 
A1  (Fig.  lb),  as  expected.  However,  close  examination 
of  the  surface  velocity  fields  reveals  that  EXP  A1  gives 
a  severely  overestimated  equatorial  current  velocity  of 
about  120  cm  s  1 ,  while  EXP  Cl  gives  a  weaker  velocity 
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of  about  70  cm  s_1.  The  observed  velocity  of  the  SEC 
is  typically  50  cm  s_1  (Reverdin  et  al.  1994),  and  it  is 
in  closer  agreement  with  EXP  Cl.  In  addition,  the  equa¬ 
torial  currents  are  not  zonal  but  show  strong  poleward 
components  in  EXP  Al  due  to  the  excessive  upwelling. 
These  unrealistic  strong  poleward  components  of  the 
equatorial  current  also  disappear  in  EXP  Cl. 

An  ideal  OMLM  should  reproduce  not  only  SST  but 
also  the  internal  structure  of  temperature  below  the  sea 
surface  in  an  OGCM.  The  most  significant  element  rep¬ 
resenting  this  is  the  mixed  layer  depth.  Figures  13-16 
show  the  distribution  of  MLD  both  from  the  observation 
(Levitus  1982)  and  the  simulations  (EXPs  Al,  Bl,  and 
Cl).  Here  MLD  is  defined  as  the  depth  of  the  temper¬ 
ature  difference  of  1°C  from  SST,  similarly  to  Rosati 


and  Miyakoda  (1988),  Ezer  (2000),  and  Sterl  and  Kat- 
tenberg  (1994).  Note  that  the  MLD  determined  in  this 
way  is  useful  to  compare  subsurface  temperature  struc¬ 
tures,  but  it  can  differ  from  a  more  physically  relevant 
MLD  based  on  the  density  criterion  in  the  area  where 
the  salinity  effects  are  important,  such  as  in  the  high- 
latitude  ocean.  Note  also  that  the  definition  of  MLD 
here  is  different  from  that  of  h  in  section  2. 

Neither  EXP  Al  nor  EXP  Bl  could  reproduce  the 
detailed  pattern  of  MLD  originating  from  the  surface 
mixing  and  circulation  of  the  ocean.  EXP  Al  produces 
a  too  shallow  mixed  layer  during  summer  owing  to  the 
absence  of  the  well-mixed  upper  layer,  whereas  EXP 
B 1  shows  very  poor  results  around  the  ACC,  as  expected 
from  the  SST  anomaly  (Fig.  7). 
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(a) Hawaii- Tahiti  Shuttle  Experiment  (Wyrtki  and  Kilonsky  1984) 


Fig.  17.  Cross  sections  of  the  temperature  (°C):  (a)  observation  [Hawaii-Tahiti  Shuttle 
Experiment  (Wyrtki  and  Kilonsky  1984)],  (b)  EXP  Al,  (c)  EXP  Bl,  and  (d)  EXP  Cl. 


Meanwhile,  EXP  Cl  shows  close  correlation  with  ob¬ 
served  data  in  the  distribution  of  MLD.  It  is  observed, 
however,  that  the  simulated  MLD  tends  to  be  slightly 
deeper  in  general.  On  the  other  hand,  the  simulated 
MLD  in  EXP  Cl  is  not  sufficiently  deep  in  the  northern 
Pacific,  which  may  be  related  to  employment  of  the 
restoring  boundary  condition  for  salinity.  Since  fresh¬ 
water  flux  plays  an  important  role  in  the  North  Pacific, 
the  restoring  boundary  condition  for  salinity  may  sup¬ 
press  convective  motion,  similarly  to  the  case  for  tem¬ 
perature.  A  noticeable  discrepancy  is  also  observed  near 
the  ACC  because  of  the  very  weak  stratification  and 
inaccurate  representation  of  the  ACC  in  this  region.  In 
particular,  an  insufficient  MLD  during  winter  in  this 
region  might  be  caused  by  excessive  diapycnal  mixing 
due  to  nonisopycnal  mixing  and  the  lack  of  ice  dynamics 
in  the  present  simulation. 

Cross  sections  of  temperature  also  illustrate  the  in¬ 
ternal  structure  of  temperature  below  the  sea  surface 
(Lig.  17).  Here  the  vertical  section  is  confined  to  the 
upper  400  m  of  the  equatorial  region  where  the  influence 
of  vertical  mixing  is  significant.  The  result  from  EXP 
Cl  shows  a  more  realistic  temperature  profile  that  re¬ 
produces  a  distinctive  thermocline  and  more  uniform 
mixed  layer,  in  agreement  with  observation  data  (Wyrtki 
and  Kilonsky  1984).  The  results  from  EXP  Al  do  not 
have  a  well-mixed  upper  layer  and  show  a  too  diffused 
thermocline.  In  EXP  Bl  temperature  decreases  too  rap¬ 
idly  below  the  EUC,  which  again  indicates  improper 
mixing  below  the  EUC,  as  in  Fig.  12. 

Finally,  Fig.  18  represents  seasonal  variations  of  the 
hemispheric-mean  heat  flux,  which  is  related  to  Fig.  9. 
It  clearly  shows  that  the  introduction  of  the  new  OMLM 
reduces  the  contribution  from  the  restoring  term  in  the 
heat  flux  substantially,  and  is  thus  expected  to  mitigate 
problems  caused  by  the  restoring  boundary  condition. 


c.  Effects  of  the  restoring  timescale  (EXPs  Cl  and 

C2) 

We  have  demonstrated  throughout  the  previous  sec¬ 
tion  that  more  realistic  features  of  the  upper  ocean,  such 
as  SST,  MLD,  and  the  equatorial  circulation,  could  be 
reproduced  by  including  the  new  OMLM.  Now  we  can 
expect  that,  if  the  OMLM  really  represents  the  vertical 
mixing  process  in  the  ocean  correctly,  it  should  repro¬ 
duce  the  observed  SST  under  a  restoring  timescale  that 
is  consistent  with  the  actual  feedback  between  the  at¬ 
mosphere  and  the  ocean.  That  is,  hopefully,  we  do  not 
have  to  use  the  unrealistically  small  restoring  timescale 
required  to  guarantee  that  the  simulated  SST  does  not 
diverse  from  the  observed  SST  in  the  previous  OGCM. 
So,  in  EXP  C2  we  attempted  to  see  how  the  results  from 
EXP  Cl  will  be  affected  when  the  restoring  timescale 
is  increased  to  a  more  realistic  value  (r  =  60  days). 

Figure  19  shows  the  global  distribution  of  SST  anom¬ 
aly  from  EXP  C2.  As  expected,  the  simulated  SST  de¬ 
viates  further  from  the  observed  SST  with  weaker  re¬ 
storing  force.  However,  the  deviation  is  largely  restrict¬ 
ed  to  amplification  of  the  existing  anomaly  of  SST  from 
EXP  Cl  (Fig.  8),  which  is  caused  basically  by  the  in¬ 
accurate  simulated  circulation.  However,  the  seasonal 
bias  of  the  hemispheric  mean  SST  still  does  not  appear. 
This  suggests  that  further  improvement  in  the  predict¬ 
ability  of  SST  can  be  obtained  under  the  present  ocean 
mixed  layer  model  if  more  accurate  circulation  can  be 
reproduced.  This  is  contrary  to  the  case  of  constant  ver¬ 
tical  mixing  (EXP  Al),  where  the  seasonal  bias  is  sig¬ 
nificantly  amplified  with  the  increased  restoring  time- 
scale  (not  shown),  and  thus  the  increased  SST  anomaly 
is  expected  regardless  of  the  circulation. 

Accordingly,  the  seasonal  cycle  of  the  hemispheric 
mean  SST  in  EXP  C2  (Fig.  20)  does  not  deviate  much 
further  compared  to  the  case  of  EXP  Cl,  although  the 
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Fig.  18.  Seasonal  variations  of  the  mean  heat  flux  (W  m  -)  (cli¬ 
matology,  EXP  Al,  EXP  Bl,  EXP  Cl):  (a)  Northern  Hemisphere, 
(b)  Southern  Hemisphere,  and  (c)  equatorial  region  (5°S-5°N). 


rms  difference  of  the  simulated  SST  from  the  observed 
SST  becomes  much  larger  (not  shown).  Furthermore, 
the  contribution  to  heat  flux  by  the  restoring  term  be¬ 
comes  even  smaller  than  the  case  of  EXP  C 1  since  the 
restoring  timescale  r  in  EXP  C2  is  twice  as  large  as  that 
in  EXP  Cl  (not  shown). 

Meanwhile,  most  other  aspects  of  the  simulation  such 
as  MLD,  the  meridional  overturing  streamfunction,  the 
surface  velocity  field,  the  volume  transport  streamfunc¬ 
tion,  and  the  cross  section  of  zonal  mean  temperature 
remain  virtually  indistinguishable  from  those  from  EXP 
Cl. 

4.  Conclusions  and  discussion 

In  this  paper  we  have  shown  that  prediction  of  the 
SST  from  an  OGCM  can  be  improved  substantially  by 
including  a  new  OMLM.  When  the  heat  flux  is  given 
by  the  form  of  the  climatological  heat  flux  plus  a  re¬ 
storing  term  toward  the  climatological  SST,  more  ac¬ 
curate  prediction  of  SST  helps  to  reduce  the  magnitude 
of  the  restoring  term,  and  thus  relieve  various  problems 
associated  with  the  restoring  boundary  condition  in  the 


OGCM.  This  implies  that  the  OGCM  with  a  new  OMLM 
reproduces  a  more  natural  ocean  along  with  more  re¬ 
alistic  heat  flux.  In  addition,  the  new  OMLM  is  found 
to  improve  various  other  aspects  of  the  OGCM,  such  as 
MLD  and  the  equatorial  circulation.  Finally,  it  is  shown 
that  the  new  OMLM  still  maintains  an  improved  level 
in  the  resultant  surface  heat  flux  and  seasonal  cycle  of 
the  hemispheric  mean  SST  under  a  longer  restoring 
timescale  of  60  days. 

Significance  of  the  present  result  can  be  found  in  the 
fact  that  it  may  provide  a  more  suitable  OGCM  for  the 
coupling  with  an  AGCM.  We  hope  that  accurate  pre¬ 
dictions  of  heat  flux  and  SST  help  remove  the  necessity 
for  flux  correction  (Meehl  1993),  and  the  reduced  ar¬ 
tificial  forcing  at  the  sea  surface  will  help  avoid  the 
spurious  oscillating  behavior  when  coupled  with  the 
AGCM  (Cai  and  Chu  1996).  Moreover,  it  is  possible 
that  the  present  OGCM  may  help  resolve  the  problem 
of  suppressed  eddy  kinetic  energy  in  high-resolution 
models  (Stammer  et  al.  1996)  since  the  artificial  re¬ 
storing  toward  the  low-resolution  climatological  SST  is 
relaxed.  We  hope  that  further  investigations  using  a 
high-resolution  OGCM  or  coupled  climate  model  can 
examine  whether  the  expected  improvement  can  be  ob¬ 
tained  in  the  near  future. 

Present  research  focuses  on  improvement  of  the 
OGCM  only  by  improving  the  vertical  mixing  process. 
Nonetheless,  there  are  many  other  parameterizations  in 
the  OGCM  that  are  important  for  more  realistic  simu¬ 
lation  of  the  ocean  (see,  e.g.,  Chassignet  and  Verron 
1998).  Certainly  more  elaborate  lateral  mixing  can  im¬ 
prove  the  simulation  results  further,  especially  in  the 
deep  convective  region  and  the  western  boundary  cur¬ 
rent  region.  On  the  other  hand,  the  present  research 
concerns  mainly  the  simulation  of  the  upper  ocean  based 
on  the  results  from  a  relatively  short  period  of  integra¬ 
tion  (24  years).  It  will  be  interesting  to  investigate  how 
the  vertical  mixing  scheme  affects  the  deep  ocean  using 
simulation  results  from  a  much  longer  period  of  inte¬ 
gration. 

Although  prediction  of  SST  was  significantly  im¬ 
proved  by  introducing  the  new  OMLM,  the  simulated 
SST  still  shows  a  slight  tendency  to  underestimate  SST 
in  the  Northern  Hemisphere  and  to  overestimate  SST 
in  the  Southern  Hemisphere  during  summer  (Fig.  9).  In 
addition,  the  simulated  MLD  tends  to  be  slightly  deeper 
than  observed.  Further  improvement  may  be  required 
in  both  the  model  and  the  forcing  data  for  more  accurate 
prediction.  The  values  of  m  in  (13)  and  z0  in  (14),  which 
determine  the  magnitude  of  eddy  diffusivity  at  the  sea 
surface,  may  need  elaboration,  utilizing  new  informa¬ 
tion  from  the  work  on  wave  breaking  and  Langmuir 
circulation.  One  also  needs  to  parameterize  the  effects 
of  high  frequency  variability  of  the  wind  stress  to  eval¬ 
uate  m  and  z0  if  monthly  mean  wind  stresses  are  used. 
More  realistic  surface  boundary  conditions  for  both  tem¬ 
perature  and  salinity  are  also  required  for  further  im¬ 
provement.  For  example,  improvement  of  the  restoring 
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Fig.  19.  Distributions  of  the  SST  anomaly  (EXP  C2)  (°C):  (a)  annual  mean,  (b)  Jan,  and  (c)  Jul. 
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Fig.  20.  As  in  Fig.  18  but  for  the  mean  SST  (°C)  (climatology, 
EXP  Cl.  EXP  C2). 

timescale  can  be  obtained  from  the  bulk  formula  for 
heat  flux  as  given  by  (4).  Finally,  it  is  important  to  notice 
that  the  value  of  a  in  (16)  and  (17),  representing  the 
effects  of  stratification,  may  depend  on  the  grid  config¬ 
uration  of  the  model,  as  mentioned  in  section  2b. 
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